Response of non-equilibrium systems with long-range initial 

correlations 



Alan Picone and Malte Henkel 

Laboratoire de Physique des Materiaux,[] Universite Henri Poincare Nancy I, 
B.P. 239, F - 54506 Vandceuvre les Nancy Cedex, France 

Abstract 

The long-time dynamics of the d-dimensional spherical model with a non-conserved order param- 
eter and quenched from an initial state with long-range correlations is studied through the exact 
calculation of the two-time autocorrelation and autoresponse functions. In the aging regime, these 
are given in terms of non-trivial universal scaling functions of both time variables. At criticality, 
five distinct types of aging are found, depending on the form of the initial correlations, while at low 
temperatures only a single type of aging exists. The autocorrelation and autoreponse exponents 
are shown to be generically different and to depend on the initial conditions. The scaling form of 
the two-time response functions agrees with a recent prediction coming from conformal invariance. 
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1 Introduction 



The slow dynamics of non-equilibrium systems displays several characteristic features which are absent 
for systems in thermodynamic equilibrium, see |TJ, [^, ^ |J for reviews. Most notable among these are the 
break-down of the fluctuation-dissipation theorem and the aging of the system as conveniently displayed 
in the two-time autocorrelations and auto-response functions. While these features were first observed 
in glassy systems, they also occur in simple ferromagnetic spin systems without disorder. The present 
paper studies such relatively simple non-equilibrium systems. 

Usually, the system is prepared in some initial state (a completely disordered initial state of effective 
infinite temperature is common) and then brought out of equilibrium even at very long times by a rapid 
quench to some final temperature T which may be below or equal to the critical temperature T c . The 
subsequent evolution then takes place at that fixed temperature T. Main observables include the two- 
time autocorrelation function C(t,s) and the auto-response function R(t,s) where s is the waiting 
time and t the observation time. A convenient way to characterize the distance of the system from 
equilibrium is through the fluctuation-dissipation ratio [|, 0, 

X(t,s)=TR(t,s)(^^J 1 (1.1) 

At equilibrium, both C and R only depend on the time difference r = t — s and X = X{j) = 1. 
Physically, this is realized if the initial quench was made to some temperature T > T c . Then for time 
scales large compared to the (finite) characteristic time scale r eq the system relaxes exponentially fast 
towards equilibrium. 

On the other hand, if either T < T c or T = T c , an infinite spin system does not reach equilibrium on 
some finite time scale but instead undergoes either phase ordering kinetics or non-equilibrium critical 
dynamics. In both cases, two-time observables such as C = C(t, s) and R = R(t, s) depend on both the 
waiting time s and the observation time t and not merely on their difference r — t — s. This breaking^ 
of time translation invariance is usally referred to as aging and will be used in this sense from now on. 
Consequently, the fluctuation-dissipation ratio X = X(t, s) ^ 1 also becomes a non-trivial function of 
both s and t. Empirically, it is well-established PJ that the aging process is associated with dynamical 
scaling, that is in the scaling regime with times 1 < s,r long enough after the initial quench one finds 

C(t, s) ~ s- b f c (t/s) , R(t, s) ~ s'^Mt/s) (1.2) 

where a, b are non-equilibrium critical exponents and fc and /r are scaling functions. Here and in the 
sequel we always have t > s. For large arguments x = t/s 3> 1, these scaling functions typically behave 
as 

f C (x)~ X -*c/z ; fR{x)r ^ x -X R / Z (13) 

where z is the dynamical exponent and Ac, A# are the autocorrelation || [| and autoresponse expo- 
nents.0 Throughout this paper, we are only interested in this late-time regime where scaling occurs. 

Usually, the long-time behaviour after a quench from a completely disordered state is studied. In 
this case, the available evidence as reviewed in || , and based on results notably from the Glauber-Ising 



model [|1(], [IT] and the spherical model with a non-conserved order parameter [[12], [13| is consistent 



2 This is associated with a breaking of the fluctuation-dissipation theorem. The time scale on which this breaking 
occurs has been studied very thouroughly in M. 

3 The values of the exponents Ac, \r (and also a, b) depend on whether T < T c or T — T c , but we shall use the same 
notation in both cases. 
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with the autocorrelation and autoreponse exponents being equal, Ac = Xr, and, for a quench exactly 
onto criticality T = T c , with the additional relation a = b. However, exceptions to this rule exist. 
Consider the 2D XY model with a fully ordered zero-temperature initial state and quenched 'upwards' 
to some temperature T < Tkt, where Tkt is the Kosterlitz-Thouless transition temperature. A recent 



calculation by Berthier, Holdsworth and Sellitto |14j] based on the spin-wave approximation has produced 
analytical results for both the autocorrelation and autoresponse scaling functions fc(x) and /r(x) and 
found the scaling "in complete agreement with the results obtained by Godreche and Luck ]73|/ for the 
spherical model at the critical point" fl4] , p. 1809]. Indeed one has a = b = rj/2, where rj = r](T) 
is the usual static (and temperature-dependent) critical exponent, but the autocorrelation exponent 
Xq/z = rj/A and the autoresponse exponent \r/z 



1 + 77/4 |14 



are different from each other. We 
therefore ask ourselves what may generically become of the relation between Ac* and Xr for more general 
initial states than fully disordered or fully ordered ones. Beyond a case study in a given model, this 
allows a test of the generic scaling and universality properties as reviewed abvoe and should provide 
useful insight. The role of the initial conditions on aging in spin-glasses has been and continues to be 
actively studied, see |L2|, [15! an d references therein. 



Questions of this sort are best addressed first in some non-trivial exactly solvable model before be- 
ing studied further through simulations. We shall therefore examine for the kinetic spherical model de- 
scribed by a non-conserved Langevin equation (see section 2 for precise definitions) the role of long-range 
correlations in the initial state which are characterized by a power law for the Fourier transformation 
C(q) of the spin-spin correlator, in the low-momentum limit \q\ — > 



C(q) 



;i.4) 



and we shall study how the above scaling forms are affected, if at all, by varying a. We can therefore 
interpolate between a disordered initial state with a = and a fully ordered state (a = —d) as considered 
[I4| . While much is already known for T < T c [plj, we shall see that for T = T c there occur several 
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new regimes for the long-time aging behaviour which depend on the value of a. Our calculations suggest 
the new scaling relation Xc = Xr + a e Q where a e g = a e g(a) is the effective value of a which actually 
describes the long-time behaviour. We shall show that this also explains the results obtained in the 2D 



XY model [141 and referred to above. 



Furthermore, introducing a new parameter into the kinetics permits an instructive test of the asser- 
tion by Godreche and Luck that the limit fluctuation-dissipation ratio 



Xoo = lim lim X(t,s) 

s^oo t— >oo 



;i.5) 



is a universal number and we shall indeed confirm this universality (see section 3). 

Finally, we recall that recently the precise form of the autoresponse function R(t, s) was derived 
from conformal invariance ||17|| , for t,s> 1 being inside the scaling regime [18|, [I! 



R(t, s) 



r 



l+a—\ R /z 



(t-s) 



-1-a 



(1.6) 



where r is a normalization constant. As we shall see in section 4, this functional form is indeed 
recovered, but the exponents a and Xr will be functions of a. 

The structure of this paper is as follows. In the next section, to give a self-contained presentation, 
we recall the main steps of the general formalism for the exact calculation of correlation and response 
functions. In section 3, we analyse these in the scaling limit where aging occurs and find the scaling 
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functions fc{ x ) an d fn{x) as functions of the initial conditions characterized by the parameter a. 
In section 4, the physical conclusions from these calculations are drawn. In the appendix, we shall 
briefly present some results on initial conditions in the ID Glauber-Ising model such that the mean 
magnetization is non- vanishing. 



2 Formalism 



We now describe the calculation of the two-time autocorrelation and autoresponse functions in the 
exactly solvable spherical model in d spatial dimensions, for arbitrary initial conditions. Our calculation 
follows closely the standard lines established several times in the past, see |L6], |2(| ^T], |22|, ^3], |24| for 
continuum field theories and 0, [F], [13], 25, 26 1 for lattice models. The effects of specific initial conditions 
on the long-time behaviour and aging will be analysed in the next section. 

The spherical model is defined in terms of real spin variables S r attached to the sites of a d- 
dimensional hypercubic lattice and subject to the constraint 



Y,Sl=Af 



(2.1) 



where M is the total number of sites, and the usual spin Hamiltonian 7i = — Yltrr') S r S r ' where the 
sum extends over nearest neighbour pairs only. The (non-conserved) dynamics is given by the stochastic 
Langevin equation 

<! ' S; J2 S s-^d + ^t))Sr + Vr(t) (2.2) 



dt 



s(r) 



where s(r) are the nearest neighbour sites of the site r, the Gaussian white noise 7] r (t) has the correlation 

(Vr{t)Vr' (0) = 2T5 ry 5(t - t') (2.3) 



and lit) is determined by satisfying the spherical constraint (|2.1|) in the mean. By a Fourier transfor- 
mation 

f(q) = E fre- iq r , /. = (27r)- d / dq f(q)e^ r (2.4) 
r Jb 

where the integral is over the first Brillouin zone B, eq. ( [2.2|) is transformed into 

dS(q,t) 



dt 



[u(q)+d(t)]S(q,t)+rj(q,t) 



(2.5) 



where in addition, together with the |qf| — ► limit 

d 



OJ 



\q) = 2^(1" cos q^^q 2 ; (fe t) rfc, t')) = 2T{2^) d 8\q + q')5(t - t') 



(2.6) 



The formal solution is |12, [T~3|1 
exp(-uj(q)t) 



S(q,t) 



S(q,Q)+ [ dt'e^'y^T^rjiq^) 
Jo 



f g(T,t) 

which forms the basis for all subsequent calculations. 



^(T,t)=exp(2 j o 



(2.7) 
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The Lagrange multiplier g(T, t) is determined from the spherical constraint and the initial conditions. 
To see this, consider the equal-time spin-spin correlation function 



Cr-r>(t) = (Sr(t)S r >(t)) 



{2.1 



where spatial translation invariance is already taken into account. Here and in the sequel, the brackets 
denote the average over the ensemble of the initial conditions and over the thermal histories, i.e. the 
realizations of the noise i] r (t). The spherical constraint (|2.1| ) implies that 



C (t) = (S r {t)S r {t)) = 1 
In Fourier space the equal-time correlator C(q,t) is obtained from 

(S(q,t)S(q',t)) = (27r) d S d (q + q')d(q,t) 

and is given by 



C(q,t) 



exp(—2uj(q)t) 



C(q, 0) + 2T f dt' e 2uj{q)t ' g(T, t') 
Jo 



(2.9) 



(2.10) 



(2.11) 



From the spherical constraint (|2.9|) we have in Fourier space that J dq (2tt) d C(q, t) = 1 and this leads 
to the following Volterra integral equation for g(T,t), see [|1^, [yj 



g(T, t) = A(t) + 2T [ dt' fit - t')g(T, t') 
Jo 

where the two auxiliary functions f(t) and A(t) are defined as follows 



(2.12) 



f(t) 



(27T) 



J dqe- 2 ^ = (e- 4 'J (4t)) d , A(t) = J dq e- 2 ^C(q, 0) (2.13) 



and where Iq is a modified Bessel function |27| . The solution of equation ( |2.12| ) is found from a Laplace 
transformation 

poo 

dtf(t)e- pt 



and is given by 



f(p) 



g{T,p) 



A(p) 



l-2T/(p) 



(2.14) 



(2.15) 



Therefore the entire evolution of the system, starting from a freely chosen initial condition, can be 
described in terms of the properties of the functions f(t) and A{t). In particular, the initial state is 
characterized exclusively by the initial correlator C(q, 0) and these data enter explicitly only into the 
function A(t). The simplest case may be considered to be an initial state without any correlations. 
Then C r (0) = 5 r ^, therefore C(q, 0) = 1 and thus A(t) = f(t). This case has been analysed in great 
detail, see |12|, 0, |26| and in particular [|13|]. 

We postpone the analysis of the effects of initial conditions given by C(q, 0) to the next section and 
now give the general expressions for the two-time correlators and the two-time response functions. The 
two-time correlation function is defined as 



C r _ r /(t, S) = (S r (t)S r '(s)) 



(2.16) 
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where s is the waiting time, t the observation time and t > s > always. In Fourier space, it is easy to 
see that 



C(q,t,s) = C(q, s)e 



-w(q)(t-s) 



i g(T,s) 
g(T,t) 



(2.17) 



and where the expression ( 2.1 1| ) for the single-time correlator has been used. Below, we shall be mainly 
be interested in the two-time autocorrelation 



C(t,s) = C (t,s) = (2n)- d ldqC(q,t,s) 



g(T } t)g(T } s) 



A[ t -±l)+2T fds' f ( ^±1 - s > )g (T, s') 



(2.18) 



which we shall analyse in the next section. The response function is obtained in the usual way [21 



12, |l3f by adding a small magnetic field term bJi = — ^2 r h r (t)S r (t) to the Hamiltonian. This leads to 
an extra term h r (t) on the right-hand side of the Langevin equation (|2.2|) . Provided that causality and 
spatial translation invariance hold, we have to first order 



(S r (t)) = / ds 2_. R r -r'(t, s)h r r(s) + . . . 



(2.19) 



which defines the (linear) response function R r (t, s). The calculation is completely standard |21], [U| 
and we merely quote the result 



R{q,t,s) 



5{S(q,t)) 



lg(T,s) 



h r =0 



5h(q, s 

Again, we shall be mainly concerned with the autoresponse function 

't-s 



g(T,t) 



(2.20) 



R(t, s) = Ro(t, s) = (27r)- d / dq R(q, t,s) = f 



g(T,t) 



(2.21) 



Eqs. (|2.11| , |2Tl8|j2.21| ) give the single-time correlation function and the two-time autocorrelation and 
autoresponse functions, respectively, for yet arbitrary initial conditions. Together with the spherical 
constraint Q2.12|) , which fixes the function g(T, t) in terms of the auxiliary functions f(t) and A(t), these 
are the main results of this section. 



3 Analysis 

It is our aim to consider the effects of long-range correlations in the initial state on the long-time and 
aging behaviour of the model. For the long-time behaviour, only the low-momentum regime should be 



relevant, which we take to be of the form [16 



C{q,Q) = Co + c a \q\ a (3.1) 

where a is a free parameter and Co, c a are normalization constants. For a > 0, the long-time behaviour 
depends to leading order only on the first term while for a < 0, we have effetively Cq = 0, to leading 
order. For the purposes of this paper, namely the study of the role of long-range correlations in the 
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initial conditions as parametrized by a, it is sufficient to study the case Cq = 0, which we shall assume 
from now on. For a < 0, this corresponds to initial correlations of the form C r (0) ~ r -d-u f Qr j ar g e 
distances r = \r\. This initial condition should lead to the following t — > oo asymptotics of the auxiliary 
function 

A(t) ~ a a t^ d+a ^ 2 , a a = c a (2n)- d Jdu e~ 2u V (3.2) 

The case without initial correlations considered previously 0, |7|, [X3j, |(| corresponds to C(q, 0) = 1 
and is recovered if we formally set a = 0. Any differences between the scaling behaviour coming 
from the initial conditions ( |3.1| ) and for those analysed previously must related to a different long-time 
behaviour of the functions A(t) and fit) (it is well known that for t large, one has fit) ~ (87rt)~ d / 2 ). 
It is the purpose of this section to carry through the technical analysis and especially to obtain the 
explicit scaling functions for the autocorrelation, autoresponse and the fluctuation-dissipation ratio. 
The physical discussion of the results will be presented in section 4. 

The single-time correlator C r (t) has been analysed in great detail by Coniglio, Ruggiero and Zanetti 
p3fl in the context of coarse-grained field theory. 

We now analyse the long-time behaviour of the two-time correlators and response functions. For 
the questions we are interested in, namely the breaking of the fluctuation-dissipation theorem and/or 
aging effects, it is enough to restrict to temperatures at or below criticality. Indeed, for T > T c , the 
system will simply relax to equilibrium within a finite time scale r eq ~ (T — T c )~ 2u , where v is the 
known equilibrium correlation length exponent, see [13| and references therein. For the same reason, 



we restrict to d > 2 throughout such that there is always a phase transition at a non-vanishing critical 
temperature T c > 0. In addition, the early stages of the evolution of the system (notably the evolution 
of the fluctuation-dissipation ratio) will depend on whether the system was initially prepared in an 
equilibrium state or not but this information requires more knowledge on the initial state than merely 
its low-momentum behaviour eq. (|3.1|). Rather, the behaviour we are interested in is described by the 



scaling or aging limit which is reached when both the waiting time s and the observation time t become 
simultaneously large, that is 1 <C s ~ t — s. As we shall be interested in the scaling properties of C(t, s) 
and R(t, s), it is useful to work with the dimensionless scaling variable 

x=->l (3.3) 

s 



Taking the scaling limit t, s — > oo with x fixed, it is convenient to rewrite the results ([2.18 , 2.21 ) of 
the previous section as follows 



C(t,s) 



s -d/2+l 



y/g(T,s)g(T,t) 



a Q s~^ ( ^±y~ 2(d+a \ [ 1 d e 9 (T,se)(x + l-29)-^ 



o 



R(t,s) * {^r d/2 s- d '\x-iy d l\l 9 ^ ) (3.4) 

where we have neglected correction terms to the leading scaling behaviour. Next, we need the leading 
behaviour of g{T,t) for t large. 
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Table 1: Values of the exponent ip = ip(pi,d) which describes the long-time behaviour of g(T c , t) accord- 
ing to eq. ( |3.5|) and the exponent (p = l + ip + a/2m the five different scaling regimes at criticality. 



Regime 


conditions 






I 


2 < d < 4 


2 < D < 4 


-1 - a/2 





II 


4 < d 


2 < D < 4 


1 - (d + a)/2 


(4-60/2 


III 


2<d<A 


4 < D 


d/2-2 


(c/ + a)/2- 1 


IV 


4<d 


4 < D a > -2 





a/2 + 1 


V 


4<d 


A < D a < -2 





a/2 + 1 



3.1 Non-equilibrium critical dynamics 

We first consider the case when T = T c , the equilibrium critical temperature. In general, we expect the 
leading asymptotics of g(T,t) to be of the form 

g{T,t)~g {a,d)t^ ; t>l (3.5) 

where the dependence of ip = ip(a, d) must be calculated and go is a constant. We determine iJj by finding 
first from eq. ( |2.15| ) the low-p behaviour of g(T,p), which in turn is given by the small-p behaviours 
of f(p) and A(p). These may be found from the integral representations (|2.14j) using well-established 
techniques [ f2"B| , |T3| . At the critical point, the results are as follows: 

J{p) ~ A 1 -A 2 p + ... ■ d>A (3.6) 

where T c = l/(2Ai) is the equilibrium critical temperature and Ak = {27i)~ d J dq (2u(q))~ k , k = 1, 2, . . . 
which exist for d > 2k. Similarly, we find 

A{p) ~ a a T (l - i(d + a)J p-^+^l 2 ■ 0<d + a<2 

A{p) ~ B 1 + 0(p,p- 1+ ^ +a )/ 2 ) ; rf + a>2 (3.7) 

where a a is given in eq. ( p.2|) and B\ = (2n)~ d fdq C(q, 0)(2u(q))^ 1 which exists for d + a >2 (if either 
d = 2, 4 or d+a = 0, 2 additional logarithmic factors will be present, which we shall discard throughout). 
From this, g~(T,p) follows and transforming back to g(T, t), the exponent ip(a, d) describing the leading 
behaviour for t large is found. We collect the results in table and identify five regimes where the 
behaviour of g(T,t) is different. It is convenient to characterize these regimes in terms of an effective 
dimension 

D = d + a + 2 (3.8) 

the meaning of which we shall discuss in section 4. From these results, we find the following scaling 
forms 



C(t,s) = (A7r)- d/2 s- d/2+1 [s- lp Mo(x)+K (x) 
R(t,s) = (ATT)- d / 2 s- d/2 (x-l)- d/2 x~^ /2 (3.9) 

where the values of the exponent ip = 1 + ip + a/2 are also listed in table |] and 



7 



a a (47r 



,d/2 



-x 



x + 1 



-(d+a)/2 



2Tx"^ /2 / uAx + 1 - 2w)- d/2 



M (x) 
K (x) 

Furthermore, the scaling of the fluctuation-dissipation ratio can be written in the form 

( r _ 1 \~d/2 T -V/2 

X(t s) = T— - - 

{h } s-*M{x)+K{x) 

where the functions M(x) and K(x) are defined by 



(3.10) 



(3.11) 



M(x) 
K(x) 



d+a ,\ , . dM (x) 
— — + i) M (x) -x- 



dx 



1 I K (x) — x 



dK (x) 
dx 



(3.12) 



and we see that the scaling of C(t, s) and consequently of X(t, s) depends on the sign of (p. For that 
reason, the regimes IV and V have to be distinguished. 

We can now list the scaling functions for both the autocorrelation and autoresponse functions 

C(t, s) = {Aix)- d ' 2 T c s~ b f c {x) , R(t, s) = (Anr^s-^fnix) (3.13) 

together with the fluctuation-dissipation ratio. The results for the exponents a and b are given in table ^| 
below. For the response function, the scaling function simply is in all five regimes 

f R {x) = {x-l)- d ' 2 x^l 2 (3.14) 

and the values of ijj can be read off from table |l|. We shall return to this result in section 4. 
The functions /c(x) and X(t, s) are listed below: for regime I, we find 
r(l - d/2)Y{-a/2) 



fc(x) 



2l+"/2 



T(l - (d + a)/2) 
r(l - d/2)r(-a/2) 



r(l-(d + a)/2) 



x 



x -(d+a/2-l)/2 +2x l/2+a/4 / &y y - l- a /2^ + X _ 2y yd/2 

JO 

d /2-a/4+l/2 x ^ OQ 15) 



The fluctuation dissipation ratio becomes in the scaling limit a function of x only, which reads for 
a ^ -2 



X = X(x) = (x-1) 



-d/2 



r(l - d/2)r(-a/2) 



r(l-(d + a)/2) 



2 a/2 



d+a _ ^ _ 1 ] x -(d+a)/2 

x+1 2 



+2 / dyy 



-l-a/2. 



x + l-2y) 



-d/2 



a d 2y — 1 



4 2 x + 1 - 2y 



r(l -d/2)r(-a/2) 



r(l-(d + a)/2) 



2 + a 



(2x 



.a/2 



x — >■ oo 



(3.16) 



together with the leading behaviour for infinitely separated timescales x — > oo. For a = —2, we find 



X = X(x) 



1 + 



x — l\ 



x+1 



d/2 



X 



X+1 



l-d/2 N 



x — > OO 



(3.17) 



S 




Figure 1: The scaling of the fluctuation-dissipation ratio X(x) as a function of x = t/s in regime I for 
d = 3 and for several values of a. In panel a, we also include the result X^ in '(x) found for regime III 
and given in eq. ( 3.19|) . In panel b we also show the function X2dxy{x) obtained in the 2D XY model 



in the spin-wave approximation [|14[ and given in eq. 

Therefore, if a 7^ —2, the limit fluctuation-dissipation ratio has the universal value = but the 
approach towards that limit does depend on the initial condition, while for a = —2, we have X^ = 1 
signalling that an equilibrium state will be reached. 

The behaviour of the scaling function X(x) is illustrated in figure [l} In the left part (figure |l|a) we 
consider initial states with a > —2 which are more disordered than the critical equilibrium state. In all 
cases, one starts from equilibrium at equal times, since X(l) = 1. If a ^ —2, the fluctuation-dissipation 
ratio decays with x increasing. However, if a approaches the border between regions I and region III, 
the scaling function X^\x) obtained for the region I goes over smoothly into the one found for region 
III. Close to that boundary, quite large values of x are needed in order to distinguish these functions. 
On the other hand, for a = —2, the long-range properties of the initial state are the same as for the 
equilibrium state. Although the system departs initially from equilibrium, since X(x) > 1, the non- 
equilibrium short-range correlations are succesively equilibrated and the system finally arrives at the 
equilibrium value of the limit fluctuation-dissipation ratio Xoo = 1. If we now consider the case a < —2 
(figure |I]b) with an initial state more ordered than the equilibrium state, the behaviour is completely 
different. Starting from X(l) = 1, the fluctuation-dissipation ratio increases with x and encounters at 
some finite value x s = x s (d, a) a singularity. For x > x s , it is negative and rapidly decays to zero with 
increasing x. 

For regime II, we find 

fc{x) = 2(4n f 2 A 2 ( 1 - *±?) x^l'-'l 2 ( """" 



2 J V 2 

X( XS s) = ' -W-^fr-i)-* 3 /x + n'^" 2 (J _ 1)/2 

1 • ' A 2 (l-(d + a)/2)(l + (d + a)(3x+l)(2x + 2)-i)\ 2 J { ' 

where the integral A 2 was defined in section 3. Therefore for any value of x, X(xs, s) = in the scaling 
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limit s — ► oo, since d > 4 here. 
For regime III, we have 



/c(s) 



4 (X - l)-rf/2+l x l-d/4 



X = X(x) 



d-2 



1 + 



x + 1 



d-2 



.r 



(3.19) 



which reproduces the scaling functions previously found [Tj| for a = 0. The entire scaling functions 
turn out to be completely independent of a. Therefore, the scaling behaviour in this region is governed 
by the effective value a e g = in the initial condition ( |3.1|) and all scaling functions are universal. In 
particular, the limit fluctuation-dissipation ratio = lim x ^ 00 X(x) = 1 — 2/d is a universal number, 
as expected. 

We have for the case of regime IV 

fc(x) = ^((x-l)- d ^-(x + l)- d ^) 



X = X(x) = 1 + 



x 



(3.20) 



which again agrees completely with the earlier results found for a = JT3[. Again, the scaling functions 
are independent of a and thus a c g = 0. In particular X^ = 1/2 is a universal number. 

Finally, in the regime V, we find 

-(d+a)/2 



fc(x) = 2(4n) d / 2 A 2 B 



v 



x + 1 



X(xs, s) 



-(4vr)- d / 2 (x - l)- d '\x + 1) fx + \V d+a ^ 2 



A 2 B v (d + a) (2x + l) 



,l+a/2 



(3.21) 



where the constant By is given by 

By- 



duu a e 



du 



u 



2u(u) 



(3.22) 



Therefore X(xs, s) = in the scaling limit, since 1 + a/2 < 0. 
From the x — * oo asymptotics 



(3.23) 



and the known fact that z — 2 Q for the spherical model at and below T c , we read off the critical 
autocorrelation and autoresponse exponents Ac and Xr and collect the results in table 0. They will be 
discussed in section 4. 



3.2 Phase ordering 

Having found the scaling behaviour at criticality, we now turn to the ordered phase, where T < T c . 
This case was considered long ago in the context of the coarse-grained 0(n)-symmetric field theory in 
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Table 2: Values of the critical autocorrelation and autoresponse exponents a, b, Ac and Xr as defined 
in eqs. ( |3.13||T23| ) in the five scaling regimes. 



Regime 


a 




b 


Ac 




I 


d/2- 


1 


d/2 - 1 


d + a/2- 1 


d- a/2-1 


II 


d/2- 


1 


1 


(d + a)/2 + l 


(d-a)/2 + l 


III 


d/2- 


1 


d/2 - 1 


3d/2-2 


3d/2 - 2 


IV 


d/2- 


1 


d/2 - 1 


d 


d 


V 


d/2- 


1 


d/2 + a/2 


d + a 


d 



the limit n — > oo |16], pi]. The relations ( |3.4| ) remain valid, but the Lagrange multiplier g(T,t) now has 
for large times the leading behaviour 



~ fl » ^-(rf+«)/2 



M 2 

cq 



(3.24) 



where M 2 q = 1 — T/T c is the squared equilibrium magnetization. In the scaling limit s — > oo, t — ► oo 



but x = t/s fixed, one has for a < 

C(t,s) 
i?(t 5 s) 
X(t, S ) 



Ax 



-(d+a)/4 



(4 



7T 



-d/2 s -d/2r _ 1 x-d/2^.(d+a)/4 



(3.25) 



4 l-(d+a)/4 T ( x + ^ 



d/2+a+l 



,-d/a+i 



{4n) d / 2 {d + a)M? ( x - l) d / 



2+1 



For a = 0, these results were already known pT| , [3~3| , |26| (for a > 0, the constant term in (|3.1| ) is 
dominating the long-time behaviour). As expected, the fluctuation-dissipation ratio X = in the 
scaling limit. The scaling functions fc(x) and fu(x) are defined as 



C{t,s) = M 2 f c {x) , R(t,s) = (4n) 



-d/2 -1-a 



fR{x) 



and we obtain 



d 

a = 1 

2 



fc(x) 



[x + 1] 
Ax 



-(d+a)/4 



f R ( X ) = {X- l)-d/2 x (d+a)/4 



(3.26) 



(3.27) 



and formally b = if we were to compare with the scaling ( |3.13|) at criticality. From these, using again 
eq. ( |3.23|) , we have the autocorrelation and autoresponse exponents 



A 



c 



~{d + a) 



Xr = -(d-a) 



(3.28) 



These results, valid for T < T c , supplement those for the critical case T = T c given in table [| For 
arbitrary a < the value of Ac was already known A long time ago, Newman and Bray [PTJ] 

studied the case of zero waiting time s = 0. From eqs ( 2.18j,272T|) and the initial condition g(T, 0) = 1, 
one has 

fit,2) - MJ-^ 



'cq 1 - 



(3.29) 



in full agreement with their results. 

4 For the 0(n)-model with n finite, Ac = for > a > a c , where a c is known to leading order in 1/n. The result 
of eq. (3.28) for Ac only holds true if a < a c . In the n — > oo limit, a c = [^6|. 
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4 Discussion 



The long-time behaviour of the kinetic spherical model depends on the spatial dimension d and the 
parameter a which charcaterizes the initial condition eq. ( |3.1Q . We have already introduced the effective 
dimension D, see eq. (|3.8|) . In terms of these, the equilibrium correlation function C eq (r) at criticality 
and the initial correlation function C- in i(r) = C r (0) scale for large distances r = \r\ as 

C eq (r) ~ r^" 2 ) , CUr) ~ r-< D - 2 > (4.1) 

(since the equilibrium critical exponent rj = for the sperical model). Therefore, we have studied 
the influence of the fluctuations of an effectively D- dimensional system onto a model defined in d 
dimensions!] We can conclude: 



1. In the scaling limit which describes the aging of the system, the results obtained for a quench to 
a temperature T <T C below the critical point and given in eqs. ( |3.27| , |3.28| ) show that the effect of 
initial long-range correlations persists for all times of the ongoing non-equilibrium behaviour, but 
without provoking any qualitative changes, as observed long ago [^T], Remarkably enough, 
however, for times 1 < t - s C s which in the scaling limit corresponds to x ~ 1 and for any 
initial condition, the two-time autocorrelations saturate at a plateau value C(x = 1) ~ M 2 q , see 
eq. ( |3.25| ). For larger differences between the waiting time s and the observation time t, C(t,s) 
decays to zero according to a power law which now again depends on the initial conditions and 
the systems shows thus a sort of memory of its initial state. 

Surprisingly, the results for a quench precisely to T = T c are qualitatively different. They are 
collected in tables [I] and and the resulting phase diagram is shown in figure |2|. Indeed, we find five 
regimes with different aging behaviours. These regimes are distinguished by the presence (when 
2 < d < 4 and/or 2 < D < 4) of strong critical fluctuations in either the equilibrium or the initial 
state, respectively, or else their absence (when d > 4 and/or D > 4) when either the equilibrium 
or the initial state are in the mean field regime. In addition, if both d > 4 and D > 4, it is of 
relevance whether the initial correlations decay faster than those in the equilibrium state or not. If 
d = D, the system is prepared in its critical equilibrium state, the fluctuation-dissipation theorem 
is valid and no aging occurs. That is the situation of equilibrium critical dynamics. Specifically, 
we find 

(a) Only if the initial correlations are in the mean field regime, i.e. D > 4, and if in addition the 
initial correlations decay faster than in equilibrium, the system may show a non-vanishing 
value of the fluctuation-dissipation ration Xoo. In these cases (regimes III/IV), the entire 
scaling functions fc{x), fn{x) and X(x) do not depend of the initial conditions at all and 
agree with the previously known results obtained for a disordered initial state of infinite 
temperature, see fl2|, |H| ^6[. One therefore has effective initial conditions such that a e s = 0. 
That finding is in full agreement with the expected [|13| , || universality of these scaling 
functions and of in particular. 

(b) A non-trivial result for the scaling of the fluctuation-dissipation ratio X(x) is found if both 
the equilibrium and the initial states are fluctuation-dominated (regime I). The limit value 
for very large separations x = t/s —>■ oo between the waiting time s and the observation 
time t is the universal value X^ = 0, provided that a ^ —2, but the approach towards 



5 Practically, long-range initial conditions as studied here might be realized by coupling the degrees of freedom of 
the model under study to those of another system at criticality. Alternatively, one might consider a spin system with 
long-range interactions J = J(r) and use this to prepare an initial state with long-range correlations. 
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Figure 2: Kinetic phase diagram at T = T c depending on the spatial dimension d and on the effective 
dimension D = d + a + 2 where a is the exponent of the initial condition eq. ( |3.1| ). The five different 
scaling regimes I, . . . , V are shown. The dashed line (a = 0) corresponds to the case of a completely 
disordered initial state, the dash-dotted line (a = —2) corresponds to an equilibrium initial state at 
T = T c and the fully ordered zero-temperature initial state corresponds to the line D = 2. 



this limit depends on the initial condition through X(x) ~ x~^ a ^ 2 , as was illustrated in 
figure |]. It is not impossible that the rather trivial value = which is more typical of a 
low-temperature phase might be a peculiarity of the spherical model. On the other hand, if 
a = — 2, then Xoo = 1 and the system evolves towards equilibrium. 

(c) If the equilibrium state is in the mean field regime and if in addition the initial correlations 
decay slower than in equilibrium, the fluctuation-dissipation ratio X vanishes in the scaling 
limit, independently of the value of x = t/s (regimes II/V). Despite of having fixed the tem- 
perature at its critical value, this behaviour is typical for a quench into the low-temperature 
ordered phase. 

It would be of interest to see whether a similar variety of different types of non-equilibrium critical 
dynamics may be established for different spin systems (especially with 2^2) and in particular, 
whether the role of the equilibrium and/or the initial state being in the mean- field regime can be 
confirmed. As a preparation for this we briefly treat in the appendix an example where the mean 
magnetization does not vanish. 

2. Turning to the values of the autocorrelation exponent Ac and the autoresponse exponent A#, we 
find that one of the following two scenarios is realized: either (i) the two exponents are equal 

Ac = Xr (4.2) 

which occurs at criticality in the regimes III and IV characterized by a c s = and where also 
Xoo ^ or else (ii) they satisfy the relation 

Ac = A fl + a cS (4.3) 
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which is realized in the entire ordered phase and for the critical regimes I, II and V, where 
a c s = a and Xoo = 0. As an extreme case, this also includes a completely ordered initial state, 
where C- m i(r) ~ cste. which corresponds to a = —d. Up to eventual logarithmic corrections, this 
case may be included by taking the limit D — > 2+. 

We point out that the results Xc/z = rj/4 and Xr/z = 1 + r//4 obtained for the 2D XY model 



14| for a fully ordered initial state (therefore a = —d = —2) in the low-temperature phase are 
consistent with eq. ( |4.3j ), since z = 2 in that model. Indeed, the behaviour of the fluctuation- 
dissipation ratio X2dxy(x) of the 2D XY model is quite analogous to the one found here for the 
spherical model in regime I at criticality. In figure |I|b we show the fluctuation-dissipation ratio of 
the 2D XY model when starting form a fully ordered initial state 



X 2 dxy(x) 



1 - 



2(x+ 1) 



(4.4) 



which is valid in the entire low-temperature phase, using the spin- wave approximation [|14]]. The 
similarity with the functions X(x) of the spherical model in regime I with a < —2 is evident. 

Therefore, we conjecture that the relationship between Ac* and Xr is in general given by eq. ( |4.3| ) 
and that only in those special cases where a e s = applies these two exponents happen to be 
equal. That was indeed the case in all models reviewed in Q]. Tests of the conjectured scaling 
relation (14. 31) in other models would be most welcome. 



Similarly, the usual anticipation, see e.g. [|J, that at criticality the two exponents a = b = 2/3 /uz, 
where (3 and v are standard equilibrium critical exponents, can be checked through the entries of 
table 0. It appears again more as a property of certain initial conditions than as a generally valid 
statement. From table 0, a — b appears to hold whenever X(t, s) does not vanish identically for 
all x in the scaling limit. In these cases the proposed relationship to the exponents (3, u, z seems 
rather to be a hyperscaling relation since it does not hold in regime IV. 

3. Having examined the asymptotic properties of the scaling functions fc{x) and Jr(x) for x large, 
we now turn to their functional form for finite values of x. Indeed, as already mentioned in the 
introduction, conformal invariance predicts that the scaling function of the autoresponse function 



should be given by [IS, IS 



f R (x) = x l+a -^ z (x - I)' 1 -* (4.5) 

That prediction can be tested by comparing with the exact result fl3.14j ) for a critical quench and 
with ( |3.27| ) for a quench into the orderd phase. Inserting the values of the exponents a and Xr 
which can be read from table |2| and from eqs. (|3.27| , |3T28| ), respectively, show perfect agreement, 



both below criticality and at criticality for all five aging regimes. 

In addition, we may also test the full space-time dependent response function. Since the dynamical 



exponent z = 2 in our model, conformal invariance predicts for z = 2 [18|, |29 



/ M r 2 \ 

R r (t, s) = R(t, s) exp (- — —) (4.6) 

where spatial translation invariance is already taken into accout, the autoresponse function R(t, s) 
is given as before by eq. ( |1.6D and M. is a non- universal and dimensionful constant. To check this, 
we calculate the full response function from eq. ( |2.20| ). In the scaling limit we are interested in, 



where both t, s as well as their difference t — s become simulataneously large, we may use the 
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q — ► limit (|2.6|) in the integral 



R r (t, s) = ( 27 r)- d ^m | dg exp [-u;(<z)(t - a ) - i<? • r] (4.7) 

and it is easy to see that the resulting gaussian integrals reproduce eq. ( |4.6|) exactly, with M. = 1/2. 

This suggests that the presence of conformal invariance in aging systems should be independent 
of spatially long-ranged correlations in the initial state. Further tests of conformal invariance in 
different aging systems with spatially long-range initial conditions are called for. 

Summarising, the study of the influence of spatially long-range correlations in the initial state 
on the long-time behaviour of two-time observables of the exactly solvable spherical model has led 
us to the recognition of several new types of non-equilibrium critical dynamics in that model. As a 
consequence, we could formulate the conjecture eq. ( |4.3|) on the relationship between the autocorrelation 
exponent Ac and the autoresponse exponent \r. Finally, the hypothesis of conformal invariance in aging 
ferromagnetic systems could be tested in a new way. 
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Appendix. Some remarks on the ID Glauber-Ising model 



Here we consider briefly some non fully disordered initial conditions in the ID kinetic Ising model with 
Glauber [|3D| dynamics. The Hamiltonian is on a chain of iV sites with periodic boundary conditions 

N 

n = -jj^wi+i ( A1 ) 

i=l 

where <7j = ±1 are the Ising spins. The dynamics may be given through the heat bath rule, which gives 
the probability of finding the spin variables <Tj(t + 1) in terms of those at time t 



P{a i {t + l) = ±l)= l - 



1 ± ta.uk ( -(<7i_i(t) + <T i+1 (t)) 



(A2) 



Here we consider the following initial conditions in terms of the spin-spin correlator C r (t) = (a r (t)a (t)) 

0.(0) = Ml + (1 - Ml) 5 r , (A3) 

where Mq = M r (0) is the initial averaged magnetization. This allows a simple case study of situations 
where the mean magnetization M r (t) does not vanish. Here and in the sequel spatial translation 
invariance is already taken into account. We study the long-time evolution of the ID model after a 
quench to zero temperature at time t = 0. 

The exact solution of the model is closely parallel to the standard lines as presented for the special 
case Mo = in flO], |Tl| . We shall therefore merely state our results. First, the mean magnetiza- 
tion M r (t) = M r (0) = Mo does remain constant for all times t. We are interested in the connected 
autocorrelation function T(t, s) and the autoresponse function R(t, s) defined as 

T(t,s) = (a r (t)a r (s))-(a r (t))(a r (s)) (A4) 
6M r (t) 



R(t,s) = T 



H r =0 



SH r (s) 

where H r is an external magnetic field at the site r. 

In ID, we find in the scaling limit t, s —>■ oo but with x = t/s > 1 fixed 



(A5) 



r(t, s) = (1 - Ml) 1 arctan yj , R(t, s) = (l - M 2 ) -U ^J-_ ( A6 ) 

and the fluctuation-dissipation ratio X(t, s) = R(t, s)(dT(t, s)/ds)~ 1 becomes a function of x only and 
reads X(t,s) = X(x) = (x + l)/2x. Up to the prefactor 1 -M r (t)M r (s) = 1 - M 2 , the results for T(t, s) 
and R(t, s) (and therefore also for X(x)) are exactly the same as found in [|K], [□]] for the autocorrelation 
and autoreponse functions at M = 0. Therefore, unless M = 1 and the system is prepared in a zero- 
temperature equilibrium state, the exponent a introduced in the text takes the value a c g = here. 
That agrees with the fact that the autocorrelation and autoresponse exponents are equal, Ac = Xr — 1. 

It would be interesting to study the effects of non fully disordered initial conditions in wider settings, 
e.g. in higher dimensions, conserved order parameters and so on. 
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